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Systems of polynomial equations over the complex or real numbers can be used to model combinatorial 
problems. In this way, a combinatorial problem is feasible (e.g. a graph is 3-colorable, hamiltonian, etc.) 
if and only if a related system of polynomial equations has a solution. In the first part of this paper, we 
construct new polynomial encodings for the problems of finding in a graph its longest cycle, the largest 
planar subgraph, the edge-chromatic number, or the largest fc-colorable subgraph. 

For an infeasible polynomial system, the (complex) Hilbert Nullstellensatz gives a certificate that the 
associated combinatorial problem is infeasible. Thus, unless P = NP , there must exist an infinite sequence 
of infeasible instances of each hard combinatorial problem for which the minimum degree of a Hilbert 
Nullstellensatz certificate of the associated polynomial system grows. 

We show that the minimum-degree of a Nullstellensatz certificate for the non-existence of a stable set of 
size greater than the stability number of the graph is the stability number of the graph. Moreover, such a 
certificate contains at least one term per stable set of G . In contrast, for non-3-colorability, we found only 
graphs with Nullstellensatz certificates of degree four. 



1. Introduction 

N. Alon pQ used the term "polynomial method" to refer to the use of non-linear polynomials for solving 
combinatorial problems. Although the polynomial method is not yet as widely used by combinatorists as, for 
instance, polyhedral or probabilistic techniques, the literature in this subject continues to grow. Prior work on 
encoding combinatorial properties included colorings [U HJ [TU1 [2H [^J [2H1 US] , stable sets [51 |23[ [55] . 
matchings [IT], and flows [2J HH1 EH]- Non-linear encodings of combinatorial problems are often compact. 
This contrasts with the exponential sizes of systems of linear inequalities that describe the convex hull 
of incidence vectors of many combinatorial structures (see [37]). In this article we present new encodings 
for other combinatorial problems, and we discuss applications of polynomial encodings to combinatorial 
optimization and to computational complexity. 

t Research supported in part by an IBM Open Collaborative Research Award and by NSF grant DMS-0608785 

t Research supported by the ISF (Israel Science Foundation) and by the fund for the promotion of research at Technion 
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Recent work demonstrates that one can derive good semidefinite programming relaxations for combina- 
torial optimization problems from the encodings of these problems as polynomial systems (see [25] and 
references therein for details). Lasserre [2D], Laurent [2T] and Parrilo [3TJ [32] studied the problem of min- 
imizing a general polynomial function f{x) over an algebraic variety having only finitely many solutions. 
Laurent proved that when the variety consists of the solutions of a zero-dimensional radical ideal / , there 
is a way to set up the optimization problem min{/(x) : x G variety (/)} as a finite sequence of semidefinite 
programs terminating with the optimal solution (see |21J. 

This immediately suggests an application of the polynomial method to combinatorial optimization prob- 
lems: Encode your problem with polynomials equations in K[a;i, . . . ,x n ] that generate a zero-dimensional 
(variety is finite) radical ideal, then generate the finite sequence of SDPs following the method in [21] . This 
highlights the importance of finding systems of polynomials for various combinatorial optimization prob- 
lems. The first half of this paper proposes new polynomial system encodings for the problems, with respect 
to an input graph, of finding a longest cycle, a largest planar subgraph, a largest fc-colorable subgraph, or a 
minimum edge coloring. In particular, we establish the following result. 

Theorem 1.1. 

1. A simple graph G with nodes 1, . . . , n has a cycle of length L if and only if the following zero-dimensional 
system of polynomial equations has a solution: 

n 

For every node i = 1, . . . , n: 



(lh-l)=0, f[{xi-a)=0, (1.2) 



jeAdj(i) 

Here Adj(i) denotes the set of nodes adjacent to node i . 
2. Let G be a simple graph with n nodes and m edges. G has a planar subgraph with K edges if and only if 
the following zero- dimensional system of equations has a solution: 
For every edge G E(G): 

z {ij} ~ z fe> = °> -K = . (1.4) 

{i,j}GE(G) 

For k = 1,2,3 , every node i G V(G) and every edge {i,j} £ E(G): 

W (x{i}k - s) = 0, JJ (y{ij} k - s) = , (1.5) 

s=l s = l 



Sk 



Yl (%{i}k - X {j}k) i i (X{i}k - V{uv}k) Yl (y{iJ}k - V{uv}k) 

1 .,jev(G) iev(G), {i,j},{u,v}eE(G) , 

\ l <3 {u,v}£E(G) ) 



1 . (1.6) 



For k = 1,2, 3 , and for every pair of a node i e V(G) and incident edge {i,j} G E(G): 

Z {ij}(y{ij}k ~ x {i}k - = . (1.7) 
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For every pair of a node i £ V(G) and edge {u, v} G E(G) that is not incident on i: 

z {uv} (V{uv}l - x {i}l - A{uv,i}l) (V{uv}2 - x {i}2 ~ A{ UOji } 2 ) (V{uv}3 ~ x {i}3 ~ A {uv,i}3) = 7 (1-8) 

z {uv} ( x {i}l - V{uv}l - A{i,uv}l) ( x {i}2 - V{uv}2 ~ A{i,uv}2) ( x {i}3 ~ V{uv}3 ~ ^{i,uv}3J = ■ (1-9) 

For every pair of edges {u,v} € E{G) (regardless of whether or not they share an endpoint): 

Z {ij} z {uv}(y{ij}l-U{uv}l~ A{ij,uv}l) (y{ij}2 — y{uv}2 — A{ij,uv}2) (y{ij}3-y{uv}3- A{ij,uv}3) = , (1.10) 

z {ij} z {uv}(y{uv}i-y{ij}i-^{uv,ij}i)(y{uv}2-y{ij}2-^{uv,ij}2)(y{uv}3-y{ij}3-^{uv,ij}3) = ■ (l.ll) 

For every pair of nodes i,j £ V(G) , (regardless of whether or not they are adjacent): 

(x {i}1 - x {j}1 - A {lj}1 ) (x {t}2 - x {j}2 - A {lj}2 ) (x {l}3 - x {j}3 - A {iJ}3 ) = , (1.12) 

( X {J}1 ~ X {i}l ~ A {j,i}l) ( X {J}2 " X { 1 }2 - A -, i}2 ) ( x {j}3 ~ X{i}3 ~ A {j,i}s) = ■ ( L13 ) 

For every Ai nc iex (e-9-, ^{ij,uv}kt A {ij,i}k > etc.) variable appearing in the above system: 

n-\-m— 1 

Yl (Aindex ~ d) = . 



d=l 



3. A graph G has a k-colorable subgraph with R 
of equations has a solution: 



(1.14) 

if and only if the following zero-dimensional system 

(1.15) 



For every vertex i G V(G): 



For every edge {i,j} € E{G): 



y% 



E vn-R = 0- 

{i,j}eE(G) 

Vij = 0, Vij {^r 1 + x t 2x 3 + ■■■ + ^f 1 ) = 



(1.16) 



(1.17) 



4- Let G be a simple graph with maximum vertex degree A . The graph G has edge- chromatic number A if 
and only if the following zero-dimensional system of polynomials has a solution: 
For every edge {i,j} € E{G): 

4 = 1. (1.18) 

For every node i G V(G): 

( \ 



J^J ( x ij x ik) 



j,k£Adj(i) 

\ j<k 



(1.19) 



/ 



where Adj(i) is the set of nodes adjacent to node i . 

[By Vizing's theorem, if the system has no solution, then G has edge- chromatic number A + 1 .] 



In the second half of the article, we look at the connection between polynomial systems and computational 
complexity. We have already mentioned that semidefinite programming is one way to approach optimization. 
It is natural to ask how big are such SDPs. For simplicity of analysis, we look at the case of feasibility instead 
of optimization. In this case, the SDPs are replaced by a large-scale linear algebra problem. We will discuss 
details in Section [3] For a hard optimization problem, say Max-Cut, we associate a system of polynomial 
equations J such that the system has a solution if and only if the problem has a feasible solution. On 
the other hand, the famous Hilbert Nullstellensatz (see [7]) states that a system of polynomial equations 
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J = {fi(x) = 0, f2{x) = 0, . . . , f r (x) = 0} with complex coefficients has no solution in C" if and only if there 
exist polynomials a\, . . . , a r € C[x±, . . . , x n ] such that 1 = ^2 &ifi ■ Thus, if the polynomial system J has 
no solution, there exists a certificate that the combinatorial optimization problem is infeasible. 

There are well-known upper bounds for the degrees of the coefficients on in the Hilbert Nullstellensatz 
certificate for general systems of polynomials, and they turn out to be sharp (see |18j). For instance, the 
following well-known example shows that the degree of ol\ is at least d m : 

fl X^, f'2 X\ , . . . , fm—1 — X rn —2 ^m— 1' f m = 1 X m _iX m 

But polynomial systems for combinatorial optimization are special. One question is how complicated are 
the degrees of Nullstellensatz certificates of infcasibility? As we will see in Section [3] unless P = NP , for 
every hard combinatorial problem, there must exist an infinite sequence of infeasible instances for which the 
minimum degree of a Nullstellensatz certificate, for the associated system of polynomials, grows arbitrarily 
large. This was first observed by L. Lovasz who proposed the problem of finding explicit graphs in A 
main contribution of this article is to exhibit such growth of degree explicitly. In the second part of the paper 
we discuss the growth of degree for the NP-complete problems stable set and 3-colorability. We establish the 
following theorem: 



Theorem 1.2. 

1. Given a graph G , let a(G) denote its stability number. A minimum- degree Nullstellensatz certificate for 
the non-existence of a stable set of size greater than u{G) has degree equal to a{G) and contains at least 
one term per stable set in G . 

2. Every Nullstellensatz certificate for non- 3-colorability of a graph has degree at least four. Moreover, in 
the case of a graph containing an odd-wheel or a clique as a subgraph, a minimum- degree Nullstellensatz 
certificate for non- 3-colorability has degree exactly four. 



The paper is organized as follows. Our encoding results for longest cycle and largest planar subgraph 
appear in Subsection |2.1| As a direct consequence, we recover a polynomial system characterization of the 
hamiltonian cycle problem. Similarly, we discuss how to express, in terms of polynomials, the decision question 
of whether a poset has dimension p . The encodings for edge-chromatic number and largest fc-colorable 
subgraph also appear in Subsection |2.1| As we mentioned earlier, colorability problems were among the first 
studied using the polynomial method; we revisit those earlier results and end Subsection |2.2| by proposing a 
notion of dual coloring derived from our algebraic set up. In Section [3] we discuss how the growth of degree 
in the Nullstellensatz occurs under the assumption P NP . We also sketch a linear algebra procedure 



we used to compute minimum-degree Nullstellensatz certificates for particular graphs. In Subsection 3.1 



we demonstrate the degree growth of Nullstellensatz certificates for the stable set problem. In contrast, in 
Subsection 3.2 we exhibit many non-3-colorable graphs where there is no growth of degree. 



2. Encodings 

In this section, we focus on how to find new polynomial encodings of some combinatorial optimization 
problems. We begin by recalling two nice results in the polynomial method that will be used later on. D. 
Bayer established a characterization of 3-colorability via a system of polynomial equations 4 . We generalize 
Bayer's result as follows: 

Lemma 2.1. The graph G is k-colorable if and only if the following zero- dimensional system of equations 

xf — 1 = 0, for every node i £ V(G), 
x k ~ x + x^~ 2 Xj H h x]~ l = 0, for every edge {i,j} € E(G) , 
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has a solution. Moreover, the number of solutions equals the number of distinct k-colorings multiplied by k\ . 

Recall that a stable set or independent set in a graph G is a subset of vertices such that no two vertices 
in the subset are adjacent. The maximum size a(G) of a stable set is called the stability number of G . We 
view the stable sets in terms of their incidence vectors. These are 0/1 vectors of length \V\ , one for every 
stable set, where a one in the i-th entry indicates that the z-th vertex is a member of the associated stable 
set. These 0/1 vectors can be fully described by a small system of quadratic equations: 

Lemma 2.2 (L. Lovasz |24j). The graph G has stability number at least k if and only if the following 
zero- dimensional system of equations 

x"l — Xi = 0, for every node i € V(G), (2-1) 
XiXj = 0, for every edge € E{G), (2-2) 

n 

$>* = *:, (2.3) 

i=i 

has a solution. 



Example 2.3. Consider the Petersen graph labeled as in Figure [T] If we wish to check whether there are 



1 




Figure 1. Petersen graph 

stable sets of size four, we take the ideal I generated by the polynomials in Eq. |2.1| |2.2| and |2.3| 

J = {x\ - Xi, X% - X 2 , X% - X 3 , x\ - Xi, x\ - X 5 , x\ - X 6 , Xy - Xj, x\ - X S , X% - Xg, x\ - X W , 
XiX & ,X2X S ,X 3 ,X 10 , X4X7, X 5 X 9 ,XiX 2 ,X 2 X 3 , X3X4, X 4 X 5 , X1X5, X 6 X7, X7X 8 , X 8 X 9 , XgXiQ, x 6 x 10 , 

xi + x 2 + x 3 + X4 + x 5 + x 6 + x 7 + x 8 + x 9 + x w - 4) . 

By construction, we know that the quotient ring R :— C[x\, . . . , xio]/I is a finite-dimensional C-vector 
space. Because the ideal I is radical, its dimension equals the number of stable sets of cardinality four in 
the Petersen graph (not taking symmetries into account). Using Grobner bases, we find that the monomials 
1, x±o, Xg, X8, xi form a vector-space basis of R and that there are no solutions with cardinality five; thus 
a(Petersen) = 4 . It is important to stress that we can recover the five different maximum-cardinality stable 
sets from the knowledge of the complex finite-dimensional vector space basis of R (see [8]). 

Next we establish similar encodings for the combinatorial problems stated in Theorem |1.1| 

2.1. Proof of Theorem 1.1 

Proof (Theorem |1.1[ Part 1). Suppose that a cycle C of length L exists in the graph G . We set = 1 
or depending on whether node i is on C or not. Next, starting the numbering at any node of C , we set 
Xi = j if node i is the j-th node of C . It is easy to check that Eqs. |l.l| and |l.2| are satisfied. 

To verify Eq. |1.3| note that since C has length L , if vertex i is the j-th node of the cycle, then one of its 
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neighbors, say k , must be the "follower", namely the (j + l)-th element of the cycle. If j < L , then the 
factor (x-i — Xk — 1) = appears in the product equation associated with the ?-th vertex, and the product is 
zero. If j = L , then the factor (x.i — x^ — (L — 1)) = appears, and the product is again 0. Since this is true 
for all vertices that are turned "on" , and for all vertices that are "off' , we have Eq. |1.3| automatically equal 
to zero, all of the equations of the polynomials vanish. 

Conversely, from a solution of the system above, we see that L variables yt are not zero; call this set C . 
We claim that the nodes i € C must form a cycle. Since yi ^ , the polynomial of Eq. |1.3| must vanish; thus 
for some j e C , 

(xi — Xj + 1) = , or (xi — Xj — (L — 1)) = . 



Note that Eq. 1.3 reduces to this form when yi ~ \ . Therefore, either vertex i is adjacent to a vertex j (with 
yj = 1) such that Xj equals the next integer value (xi + 1 = Xj), or xi — L — Xj — 1 (again, with yj = 1). 
In the second case, since Xi and Xj are integers between 1 and L , this forces Xi = L and xj = 1 . By the 
pigeonhole principle, this implies that all integer values from 1 to L must be assigned to some node in C 
starting at vertex 1 and ending at L (which is adjacent to the node receiving 1). □ 



We have the following corollary. 



Corollary 2.4. A graph G has a hamiltonian cycle if and only if the following zero-dimensional system of 
2n equations has a solution. For every node i G V{G) , we have two equations: 

n 

W_{ x i — s ) = , and Y\ ( x i ~~ x j + ~ x j — (n — 1)) = ■ 

s=l jeAd](i) 

The number of hamiltonian cycles in the graph G equals the number of solutions of the system divided by 
2n . 



Proof. Clearly when L = n we can just fix all yi to 1 , thus many of the equations simplify or become 
obsolete. We only have to check the last statement on the number of hamiltonian cycles. For that, we remark 
that no solution appears with multiplicity because the ideal is radical. That the ideal is radical is implied by 
the fact that every variable appears as the only variable in a unique square-free polynomial (see page 246 of 
[T^]). Finally, note for every cycle there are n ways to choose the initial node to be labeled as 1 , and then 
two possible directions to continue the labeling. □ 



Note that similar results can be established for the directed graph version, thus one can consider paths 
or cycles with orientation. Also note that, we can use the polynomials systems above to investigate the 
distribution of cycle lengths in a graph (similarly for path lengths and cut sizes). This topic has several 
outstanding questions. For example, a still unresolved question of Erdos and Gyarfas [3S] asks: If G is a 
graph with minimum-degree three, is it true that G always has a cycle having length that is a power of two? 
Define the cycle-length polynomial as the square-free univariate polynomial whose roots are the possible 
cycle lengths of a graph (same can be done for cuts). Considering L as a variable, the reduced lexicographic 
Grobner basis (with L the last variable) computation provides us with a unique univariate polynomial on L 
that is divisible by the cycle-length polynomial of G . 

Now we proceed to the proof of part 2 of Theorem |1.1| For this we recall Schnyder's characterization of 
planarity in terms of the dimension of a poset |33j : For an n-element poset P , a linear extension is an order 
preserving bijection a : P — > {1, 2, . . . , n} . The poset dimension of P is the smallest integer t for which there 
exists a family of t linear extensions <Ti, . . . ,o~t of P such that x < y in P if and only if cr^x) < <7;(y) for all 
o~i . The incidence poset P{G) of a graph G with node set V and edge set E is the partially ordered set of 
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height two on the union of nodes and edges, where we say x < y if x is a node and y is an edge, and y is 
incident to x . 

Lemma 2.5 (Schnyder's theorem |33j). A graph G is planar if and only if the poset dimension of P(G) 
is no more than three. 

Thus our first step is to encode the linear extensions and the poset dimension of a poset P in terms of 
polynomial equations. The idea is similar to our characterization of cycles via permutations. 



Lemma 2.6. The poset P = (E, >) has poset dimension at most p if and only if the following system of 
equations has a solution: 
For k — 1 , . . . , p : 

H( Xi (k)-s)=0, for every i £ {1, ... ,\E\}, and 

II Xi(k)-Xj{k) 1=1. (2.4) 



s=l 



{i,i}e{i,...,|E|}, 
i<j 



For k = \, . . . ,p , and every ordered pair of comparable elements > e 3 in P : 

Xi{k) — Xj{k) — Aij(k) = 0. 

For every ordered pair of incomparable elements of P (i.e., e, ^ ej and ej j£ e.i) : 
p p 

H (x t (k) - Xj (k) - Ay(fc)) = , .[ ( Xj (k) - Xi (k) - A^k)) = , 



k=l 



k=l 



For k = 1, . . . ,p , and for every pair {i, j} £ {1, . . . , \E\}: 



\E\-1 \E\-1 

JJ (A ij -(A)-d)=Q 1 Y[(A ]t (k)-d) = 

d=l d=l 



(2.5) 



(2.6) 



(2.7) 



Proof. With Eqs. 2.4 and 2.5 we assign distinct numbers 1 through \E\ to the poset elements, such that 



the properties of a linear extension are satisfied. Eqs. |2.4| and |2.5| are repeated p times, so p linear extensions 
are created. If the intersection of these extensions is indeed equal to the original poset P , then for every 
incomparable pair of elements in P , at least one of the p linear extensions must detect the incomparability. 
But this is indeed the case for Eq. |2.6[ which says that for the l-th linear extension, the values assigned to 
the incomparable pair ei,ej do not satisfy Xi(l) < Xj(l) , but instead satisfy Xj(l) > Xi(l) . □ 



Proof (Theorem |l.l| Part 2). We simply apply the above lemma to the particular pairs of order relations 
of the incidence poset of the graph. Note that in the formulation we added variables z^j} that have the effect 
of turning on or off an edge of the input graph. □ 

Example 2.7 (Posets and Planar Graphs). 



Proof (Theorem |1.1[ Part 3). Using Lemma 2.1 we can finish the proof of Part 3. For a fc-colorable 
subgraph H of size R , we set = 1 if edge {i, j} £ E(H) or yij = otherwise. By Lemma 2.1 the resulting 
subsystem of equations has a solution. Conversely from a solution, the subgraph H in question is read off 
from those y^ ^ . Solvability implies that H is fc-colorable. □ 
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Figure 2. Via Schnyder's theorem, the square is planar since P(square) has dimension at most three. 



Before we prove Theorem |1.1| Part 4, we recall that the edge- chromatic number of a graph is the minimum 
number of colors necessary to color every edge of a graph such that no two edges of the same color are incident 
on the same vertex. 



Proof (Theorem Part 4). If the system of equations has a solution, then Eq. 1.18 insures that 
all variables Xij are assigned A roots of unity. Eq. 1.19 insures that no node is incident on two edges of 
the same color. Since the graph contains a vertex of degree A , the graph cannot have an edge-chromatic 
number less than A , and since the graph is edge-A-colorable, this implies that the graph has edge-chromatic 
number exactly A . Conversely, if the graph has an edge-A-coloring, simply map the coloring to the A roots 
of unity and all equations are satisfied. Since Vizing's classic result shows that any graph with maximum 
vertex degree A can be edge-colored with at most A + 1 colors, if there is no solution, then the graph must 
have an edge-chromatic number of A + 1 . □ 



2.2. Normal forms and Dual colorings 

In [5] Alon and Tarsi show another polynomial encoding of fc-colorability. Here we consider one curious 
consequence of the polynomial method for graph colorings when we use an algebraic encoding similar to that 
of P]. By taking a closer look at the normal form of the polynomials involved, we can derive a notion of dual 
coloring, which has the nice property that a graph is dually d-colorable if and only if it is <i-colorable. This 
gives rise to an appealing new graph invariant: the simultaneous chromatic number o~(G) , defined to be the 
infimal d such that G has a <i-labeling that is simultaneously a coloring and a dual coloring. 

Fix a graph G — (V,E) with V := {l,...,n} and E C ( 2 ) , fix a positive integer d , and let D := 
{0, 1 . . . , d — 1} . Let a := exp(^p) <E C be the primitive complex d-th root of unity, so that a , . . . , a^ 1 
are distinct and a d = 1 . For a d-labeling c : V — ► D of the vertices of G , let 

e( C ):=n{K (i) -" C(?) ) : * < J, ihj] e E } ■ 

Clearly, c is a proper d-coloring of G if and only if e(c) ^ . 

With every orientation O = (V, A) of G (where A denotes the set of "arrows" or directed edges) associate 
a sign sign = ±1 defined by the parity of the number \{(i, j) € A : i > j}\ of flips of O from the standard 
orientation (where every directed edge (i, j) has i < j), and an out-degree vector 5° :— (5°, . . . , <5°) with Sf 
the out-degree of vertex i in O . For a non-negative integer k let [k] £ D be the representative of k modulo 
d , and for a vector 5 = (5 X , ... ,S n ) eV n let [S] = ([Si], [S n ]) e D v . For a labeling c* : V — ► D of the 
vertices of G let 

e *(c*) := { sig n ° : O orientation of G with [S°] = c* } . 
Call c* a dual d-coloring of G if e*(c*) =/= . 



Theorem 2.8. A graph has a d-coloring, namely c € D v with e(c) ^0 (so is d-colorable) if and only if it 
has a dual d-coloring, namely c* € D v with e*(c*) ^0 (so is dually d-colorable). 



Proof. Let G be a graph on n vertices. Consider the following radical zero-dimensional ideal / in C[xi , . . . , x n ] 
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and its variety variety (7) in C n : 

7 := (xf-l,...,xt-l) , variety(7) := { a c := (a c(1) , . . . , a c(n) ) G C n : c G 7J y } . 

It is easy to see that the set {xf — 1, . . . , x„ — 1} is a universal Grobner basis (see [3] and references therein). 
Thus, the (congruence classes of) monomials x c , c* G 7) 17 (where 2 C := FJ™ = l x\ ), which are those 
monomials not divisible by any xf , form a vector space basis for the quotient C[xi, . . . ,x n ]/I . Therefore, 
every polynomial f = ^2a$-x s has a unique normal form [/] with respect to this basis, namely the polynomial 
that lies in the vector space spanned by the monomials x c , c* G D v , and satisfies / — [/] G 7 . It is not 
very hard to show that this normal form is given by [/] = ^2 as ■ x^ . 
Now consider the graph polynomial of G , 

fa ■= J] { (** - ^i) : * < 3, {hi} e 75 } . 

The labeling c € 7> y is a d-coloring of G if and only if e(c) = fc(a c ) ^ . Thus, G is not d-colorable if and 
only if fa vanishes on every a c G variety(7), which holds if and only if / € 7, since 7 is radical. It follows 
that G is d-colorable if and only if the representative of fa is not zero. Since /g = S s ig n ° ' 2:5 i with the 
sum extending over the 2' E ' orientations O of G , we obtain 

[ /G ] = ^sign .^°]= £ e*(c*)-* c * . 

Therefore [/ G ] ^ and G is d-colorable if and only if there is a c* G 7> y with e*(c*) ^ . □ 

Example 2.9. Consider the graph G = (V, E) having V" = {1, 2, 3,4} and E = {12, 13, 23, 24, 34} , and 
let d — 3 . The normal form of the graph polynomial can be shown to be 

[/g] = X1X2X3 ~~ X1X2X4 + X1X2X4 — X1X2X3 + x\x\xa — X1X3X4 + X\X2 — X\X2X%X4 + X1X3X4 

2 22.2 222. 22 

— X1X3 + X1X2X3X4 — X1X2X4 + X :i — X3X4 + X2X3X4 — X 2 + X2X4 — £2X32:4 . 

Note that in general, the number of monomials appearing in the expansion of /g can be as much as the 
number of orientations 2' E '; but usually it will be smaller due to cancellations that occur. Moreover, there will 
usually be further cancellations when moving to the normal form, so typically [/g] will have fewer monomials. 
In our example, out of the 2' E ' = 2 5 = 32 monomials corresponding to the orientations, in the expansion of 
fc only 20 appear, and in the normal form [fa] only 18 appear due to the additional cancellation: 

— [0:1X30:4] + [212:22:4] = —2:124 + X1X4 = . 

Note that the graph G in this example has only six 3-colorings (which are in fact the same up to relabeling 
of the colors), but as many as 18 dual 3-colorings c* corresponding to monomials x c appearing in [fc] ■ For 
instance, consider the labeling c*(l) = c*(2) = c*(4) = 0, c*(3) = 2: the only orientation O that satisfies 
[<$P] = c*(j) for all j is one with edges oriented as 21, 23, 24, 31, 34 , having sign = 1 and out-degrees 5° = 

(5° = , #p = 2 and 8° — 3 , contributing to [f G ] the non-zero term e*(c*) nt=i x j ^ = ^ ' ^i^^I^ = x i ■ 
Thus, c* is a dual 3-coloring (but, since c*(l) = c*(2) , it is neither a usual 3-coloring nor a simultaneous 
3-coloring — see below). 

Note that in this example, and seemingly often, there are many more dual colorings than colorings; this 
suggests a randomized heuristic to find a dual d-coloring for verifying d-colorability. 

A particularly appealing notion that arises is the following: call a vertex labeling s : V — > D a simultaneous 
d-coloring of a graph G if it is simultaneously a d-coloring and a dual d-coloring of G . The simultaneous 
chromatic number cr(G) is then the minimum d such that G has a simultaneous d-coloring. This is a strong 
notion that may prove useful for inductive arguments, perhaps in the study of the 4-color problem of planar 
graphs, and which provides an upper bound on the usual chromatic number x(G) . First note that, like the 
usual chromatic number, it can be bounded in terms of the maximum degree A(G) as follows. 
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Theorem 2.10. The simultaneous chromatic number of any graph G satisfies c(G) < A(G) + 1 . Moreover, 
for any G and d > A(G) + 1 , there is an acyclic orientation O whose out-degree vector 8 = (8®, . . . , <S„) 
provides a simultaneous d-coloring s defined by s(i) :— Sf for every vertex i . 

Proof. We prove the second (stronger) claim, by induction on the number n of vertices. For n — 1, this is 
trivially true. Suppose n > 1 , and let d :— A(G) + 1 . Pick any vertex i of maximum degree A(G) , and 
let G' be the graph obtained from G by removing vertex i and all edges incident on i . Let O' be an acyclic 
orientation of G' and s' the corresponding simultaneous d-coloring of G' guaranteed to exist by induction. 
Extend O' to an orientation O of G by orienting all edges incident on i away from i , and extend s to the 
corresponding vertex labeling of G by setting s(i) := 8f = d — 1 . Then O is acyclic, and therefore O is the 
unique orientation of G with out-degree vector 8 . Thus, 

e*(s) = ^{sign e : 6 orientation of G with [S e ] = s = 5° } =±1^0, 

and therefore s is a dual d-coloring of G . Moreover, if j is any neighbor of i in G , then the degree of j in 
G' is at most d — 2 , and therefore its label s'(j) = 8° (j) < d — 2 , and hence s(j) = s'(j) ^ d — 1 = s(i) . 
Therefore, s is also a d-coloring of G , completing the induction. □ 

Example 2.11 (simultaneous 4-coloring of the Petersen graph). 

According to Figure |] 8° = (2,1,0,2,0,3,1,2,3,1) . By inspection of Figure |] s(i) := 8? does indeed 
describe a valid 4-coloring of the Petersen graph. 




Figure 3. Left: A vertex labeling. Right: An acyclic orientation labeled with out-degrees. 



There are many fascinating new combinatorial and computational problems related to this new graph 
invariant, the behavior of which is quite different from that of the usual chromatic number. For instance, the 
direct analog of Brooks' theorem, which states that every connected graph with maximum degree A that is 
neither complete nor an odd cycle is A-colorable, fails: It is not hard to verify that the simultaneous chromatic 
number of the cycle G n is 2 if and only if n is a multiple of 4; thus, the hexagon satisfies a(Ce) = 3 > A(Ge) . 
Which are the simultaneous chromatic Brooks graphs, i.e. those with a(G) = A(G) ? What is the complexity 
of deciding if a graph is simultaneously <i-colorable? Which graphs are simultaneously c?-colorable for small 
d ? For d = 2 , the complete answer was given by L. Lovasz |5S] during a discussion at the Oberwolfach 
Mathematical Institute: 

Theorem 2.12. (Lovasz) A connected bipartite graph G = (A,B,E) has simultaneous chromatic number 
cr(G) = 2 if and only if at least one of \A\ and \B\ has the same parity as \E\ . 

3. Nullstellensatz Degree Growth in Combinatorics 

The Hilbert Nullstellensatz states that a system of polynomial equations {f\ (x) — 0, f2(x) — 0, . . . , f r (x) = 0} C 
C[xi, . . . , x n ] has no solution in C" if and only if there exist polynomials a±, . . . , a r £ C[x%, . . . , x n ] such 
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that 1 = a ifi ( see 0)- The purpose of this section is to investigate the degree growth of the coefficients 
a.i . In particular, systems of polynomials coming from combinatorial optimization. 

In our investigations, we will often need to find explicit Nullstellensatz certificates for specific graphs. This 
can be done via linear algebra. First, given a system of polynomial equations, fix a tentative degree for the 
coefficient polynomials at in the Nullstellensatz certificates. This yields a linear system of equations whose 
variables are the coefficients of the monomials of the polynomials a±, . . . , a r . Then, solve this linear system. 
If the system has a solution, we have found a Nullstellensatz certificate. Otherwise, try a higher degree for the 
polynomials on . For the Nullstellensatz certificates, the degrees of the polynomials on cannot be more than 
known bounds (see e.g., |18j and references therein), thus this is a finite (but potentially long) procedure to 
decide whether a system of polynomials is feasible or not. In practice, sometimes low degrees suffice to find 
a certificate. 



Example 3.1. Suppose we wish to test K 4 for 3-colorability, and we assume that the on in the Nullstellen- 
satz certificate have degree 1. After encoding with the system of polynomial equations, we "conjecture" 
that there exists a Nullstellensatz certificate of the following form 

1 = (cixi + c 2 x 2 + C3X3 + c 4 x 4 + c 5 )(xl - 1) + (c 6 xi + c 7 x 2 + c$x 3 + C9X4 + cio)(x2 - 1) 



+ (cuxi + 

+ (C21X1 + 
+ (C31»l + 

+ (cuxi + 
When we multiply out this 

1 = C\X l 



+ c 15 )(x 3 3 - 1) + (a 6 xi + ■■■ + c w ){x\ - 1) 

+ c 25 )(xj + xix 2 + x%) + (c 26 xi H h c 30 ){x\ + x 4 x 3 + x%) 

+ c 35 )(xl + x x x 4 + x\) + (c 36 x 4 H h c 40 )(xj + x 2 x 3 + x 2 3 ) 

+ c 45 )(xl + x 2 Xi + x\) + (c i6 xi H h c 50 )(a;3 + a; 3 a;4 + x\). 

certificate, we group together like powers of x\ , x 2 , x 3 , X4 as follows: 



+ C13Z3 



c s xlx 3 



+ (C35 + C45 + Cbo)xl H 1- (-C15 



• ■ + (c 2 2 + C 2 1 + C 27 + C 32 )x\x 2 

C20 - c 5 - C10). 



Because the Nullstellensatz certificate is identically 1 , this identity gives rise to the following system of linear 
equations: = cj.,0 = c 13 , = c 8 , = c 22 + c 21 + c 27 + c 32 , = c 35 + c 45 + c 50 , . . . , 1 = -c 15 - c 2Q - c 5 - c w . 
In other words, we have a large-scale sparse system of linear equations that consists only of Is and —Is. 
We implemented an exact- arithmetic linear system solver. In this example, it turns out that degree 1 is not 
sufficient for generating a Nullstellensatz certificate that is, this linear system has no solution. Ultimately, 
we discovered that degree four is required, and we were able to produce the following certificate: 



I 3 
{-X 1 



1) + f g^4 ^ gXlx 2 - g x l x 3 - g^4^1 + gXlx2Xl + -X 2 i X 3 X 1 j (x\ + X 2 X 4 + X%) 



1 4 2 3 
~X 4 + -X 4 X 2 



1 3 
-X 4 X\ 



9' 



X 4 X 2 X\ J (x 2 + X3X2 + X 3 ) + ( -X 4 + -X 4 X 2 + ~X 4 Xi + -X 4 X 2 X± J (x 4 + X3X4. + X 3 
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— -x 4 + x 4 x± — x 4 x 1 + x 1 (x 4 + x\x 4 + x 4 ) + -x 4 x 2 (x 2 + XiX 2 + x 4 ) + 



1 4 

~x 4 



X 4 X 2 (x 3 + X1X3 + x 4 ). 



Lemma 3.2. If P =/= NP , then there must exist an infinite family of graphs whose minimum- degree non- 
3-colorability Nullstellensatz certificates have unbounded growth with respect to the number of vertices and 
edges in the graph. 



Proof. Our proof is by contradiction with the hypothesis P ^ NP . Consider a non-3-colorable graph that 
has been encoded as the system of polynomial equations (xf — 1) = for i S V(G) , and (x\ + XiXj + xj) = 
for {i,j} G E(G) . Assume that every minimum-degree non-3-colorability Nullstellensatz certificate has 
deg(ai) < d for some constant d . We will show that P = NP by providing a polynomial-time algorithm 
for solving Graph-3-Coloring: (1) Given a graph G , encode it as the above system of polynomial equations, 
(2) Construct and solve the associated linear system for monomials of degree < d , (3) If the system has a 
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solution, a NuUstellensatz certificate exists, and the graph is non-3-colorable: Return no, (4) If the system 
does not have a solution, there does not exist a NuUstellensatz certificate, and the graph is 3-colorable: 
Return yes. 

Now we analyze the running time of this algorithm. In Step 1, our encoding has one polynomial equation 
per vertex and one polynomial equation per edge. Since there are 0(n 2 ) edges in a graph, our polynomial 
system has n+n 2 = 0(n 2 ) equations. Since every equation only contains coefficients ±1 and is of degree three 
or less, encoding the graph as the above system of polynomial equations clearly runs in polynomial-time. 

For Step 2, we note that by Corollary 3.2b of [34 , if a system of linear equations Ax — b has a solution, 
then it has a solution polynomially-bounded by the bit-sizes of the matrix A and the vector b (see [53] for a 
definition of bit-size). In this case, the vector b contains only zeros and ones. To calculate the bit-size of A , 
we recall our assumption that, for every a.i , deg(ai) < d for some constant d . Therefore, an upper bound 
on the number of terms in each a% is the total number of monomials in n variables of degree less than or 
equal to d . Therefore, the number of terms in each a% is 

) + (" ^ ~ 2 ) + ' ' ' + [lZ I) = 0(n d ) + O^ 1 ) + ■■■ + 0(1) = 0(n d ). 

Since there are 0(n 2 ) equations, there are at most 0(n d+2 ) unknowns in the linear system, and thus, 
0(n d+2 ) columns in A . Since the vertex equations (xf — 1) = have two terms, and the edge equations 
(x 2 + XiXj + x 2 ) — have 3 terms, there are 0(n d+2 ) terms in the expanded NuUstellensatz certificate, and 
0(n d+2 ) rows in A . Since entries in A are 0,±1 , the matrix A contains only entries of bit-size at most 2. 
Therefore, the bit-sizes of both A and b are polynomially-bounded in n , and by Theorem 3.3 of |34j . the 
linear system can be solved in polynomial-time. 

Therefore, we have demonstrated a polynomial-time algorithm for solving Graph-3-Coloring, and since 
Graph-3-Coloring is NP-Complete ([H]), this implies P = NP , which contradicts our hypothesis. Therefore, 
deg(a^) ^ d for any constant d . □ 



Thus, in the linear algebra approach to finding a minimum-degree NuUstellensatz certificate, the existence 
of a universal constant bounding the degree is impossible under a well-known hypothesis of complexity theory. 
Clearly, a similar result can be obtained for other encodings (for a generalized statement see [55]). Note that 
the linear algebra method does not rely on any property that is unique to a particular combinatorial or NP- 
Complete problem; the only assumption is that the problem can be represented as a system of polynomial 
equations. We will use it to find NuUstellensatz certificates of non-3-colorability and sizes of stable sets of 
graphs. 

We remark that this linear algebra method finds not only a NuUstellensatz certificate (if it exists), but 
it finds one of minimum possible degree. With our implementation, we ran several experiments. We quickly 
found out that the systems of linear equations are numerically unstable, thus it is best to use exact arithmetic 
to solve them. The systems of linear equations arc also quite large in practice, as the bound on the degree 
of the polynomial coefficients grows. Thus we need ways to reduce the number of unknowns. 

We will not discuss here ad hoc methods that depend on the particular polynomial system at hand (see [25] 
for methods specific to 3-colorability) , but one rather useful general trick is to randomly eliminate variables 
in the above procedure. Instead of allowing all monomials of degree < d to appear in the construction of the 
linear system of equations, we can randomly set unknowns in the linear system of equations to — e.g., set 
each variable to with probability p , independently, to get a smaller system. 

In Figure [4] we see the results of the probabilistic search for smaller NuUstellensatz certificates. On the 
cc-axis is the probability p of keeping an unknown in the linear system. Thus, if p = 0.1 , 90% of the time we 
set the unknown to 0, and only 10% of the time, we keep it in the system. For the cliques and odd wheels, we 
know that there is always a certificate of degree four. For every probability 0.1, 0.2, . . . , 1 we performed 100 
searches for a degree four certificate. For both the cliques and the odd wheels, at p = 0.1 and p — 0.2 , we 
almost never found a certificate. But for p = 0.4 , we found certificates 95% of the time. In practice this idea 
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Figure 4. Probability tests on cliques and odd wheels. 



is useful. We can reduce the number of variables in the linear system by 60%, and still find a Nullstcllcnsatz 



certificate 90% of the time. We will report the results of specific computations in Section 3.2 



3.1. The Nullstellensatz and Stable Sets of Graphs 

Lovasz [24 stated the challenge of finding an explicit family of graphs with growth in the degree of Nullstel- 
lensatz certificates. Here we solve his challenge. Our main result is stated in Theorem 1.2: For every graph 
G , there exists a Nullstellensatz certificate of degree a(G) , the stability number of G , certifying that G 
has no stable set of size > a(G) . Moreover, this is the minimum possible degree for all graphs. In what 
follows, for any graph G with stability number a(G) and integer r > 1 , we associate the infeasible system 
of polynomial equations J(G,r) : 



-(a(G) + r)+]T; 







2 = 1 



x\ — Xi — 0, for every node i € V(G) , 
XiXj = 0, for every edge {i,j} £ E(G) 

Thus, the Nullstellensatz certificate will have the general form: 



l = A(-(a(G) + r)+J2' 



i=l 



ieV(G) {i,j}£E(G) 



(3.1) 

(3.2) 
(3.3) 



(3.4) 



In the rest of this section, we will refer to the coefficient polynomials using these particular letters (i.e. A, Qi , 
etc). 



Definition 3.3. A Nullstellensatz certificate (Eq. 3.4 ) has degree d ifmax{deg(A), deg(Qi), deg(Qij)} = d 
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Lemma 3.4. For any graph G and a Nullstellensatz certificate 

l = A(-(a{G)+r) + J2 x i)+ E Qi( x t ~ x i) + E Qv( x i x i) ( 3 - 5 ) 

B 

certifying that G has no stable set of size (a(G)+r) (with r > I), we can construct a "reduced" Nullstellensatz 
certificate 

l = A'(-(a(G)+r) + J2 x i)+ E + E Q'M-'^- 

1. The coefficient A' multiplying — (a(G) + r) 4- X)»=i x i ^ as on ty square-free monomials supported on stable 
sets of G , and thus deg(^4') < a(G) . 

2. max.{deg(A),deg(Qi),deg(Qij)} = max{deg(A / ),deg(Q' i ) 1 deg(Q / i: j)} . Thus, if the original Nullstellensatz 
certificate has minimum- degree, the "reduced" certificate also has minimum-degree. 



Proof. Let / be the ideal generated by x\ — Xi (for every node i € V(G)), and XiXj ( for every edge {%, j} S 
E(G)). We apply reductions modulo / to Eq. 3.5 If a non-square-free monomial appears in polynomial A , say 



H «2 



xf fc with at least one a,- > 1 , then we can subtract the polynomial xf 1 ^" 2 • • • , a;- 3 x? k B(xj 



.) from AB and simultaneously add it to ^2Q s {x 2 s — x s ) . Thus, eventually we obtain a new certificate 
that has only square-free monomials in A' . Furthermore, if Q' s has new monomials, they are of degree less 
than or equal to what was originally in A . 

Similarly, if x ix Xi 2 ■ ■ ■ Xi k appears in A , but x^x^ ■ ■ ■ Xi k contains an edge {i,j} € E{G) (if XiXj divides 
Xi^i 2 ■ ■ ■ Xi k ), then we can again subtract B{xi 1 Xi 2 ■ ■ ■ x ik /xiXj){xiXj) from AB , and, at the same time, add 
it to j}eE(G) Qij x i x j ■ Furthermore, the degree is maintained, and we have reached the form we claim 
exists for A' . □ 



We now show that, for every graph, there exists an explicit Nullstellensatz certificate of degree a(G) . 

Theorem 3.5. Given a graph G , there exists a Nullstellensatz certificate of degree a(G) certifying the 
non-existence of stable sets of size greater than ot{G) . 



Proof. The proof is an algorithm to construct the explicit Nullstellensatz certificate for the non-existence 
of a stable set of size a(G) + r, with r > 1 . First, let us establish some notation: Let S(i, G) be the set 
of all stable sets of size i in G . We index nodes in the graph by integers, thus stable sets in G are subsets 
of integers. When we refer to a monomial Xd x Xd 2 • • ■ Xd t as a "stable set" , we mean {d\, . . . , di} € S(i, G) . 
We use "hat" notation to remove a variable from a monomial, meaning Xd 1 Xd 2 ' ' ' Xd i _ 1 x d ~ i Xd i+1 ■ ■ 'Xd k — 
x dl x d2 ■ ■ ■ • • ■x,i k . Finally, let 

&>{i,G):= Yl x dl x d2 ■ ■ ■ x dz , and 0>(O, G) := 1 . 

{d 1 ,d 2 ,...,d t }eS(i,G) 

When £P(i, G) is ordered lexicographically, we denote G) as the j-th term. We also define the constants 

G G — lG ^~ x and G G — 1 

L^ a{G) + r-z ' C ° - a{G) + r ' 

Our main claim is that we can construct explicit coefficients A, Qi, Qij of degree less than or equal to a(G) 
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ALGORITHM (Nullstellensatz certificate construction) 

INPUT: A graph G{V,E), associated polynomials @>{i,G) , a(G),r 

OUTPUT: polynomials (A, Qi,Qij) such that AB + C = 1 is true 



A^-J2^C?^(i,G) 

1 Qj«-0,fori = l...|y| 

2 Q^^O , for {i,j}£E(G) 

3 for i <- to a(G) 

4 for j <— 1 to # of monomials in G) 

5 for k <- 1 to |V| 

6 Let ^(i, G) = Xd-tXdz • • -Xd t 

7 if 3 3 j(i,G)xk is a square-free stable set 

8 (rule 1) Q k <- Q fc + Ci +l Xd x Xd 2 ■ ■ ■ x di 

9 else if 3?j(i,G)xk is a square- free non-stable set 

10 Choose an edge I, 1 < I < i, so that (fc, d;) is an edge of G 

12 (rule 2) Q kdl <— Q kdl + Cfx^x^ ■ ■ ■ x dl _ 1 x2 l x dl+1 ■ ■ ■ x di 

13 end if 

14 end for 

15 end for 



16 end for 

17 return A,Qi,Qij 

Figure 5. Nullstellensatz certificate construction 
such that the following identity is satisfied 

(n \ n 

- (a(G) +r) + J2 x i)+ QijXiXj + ^Qi{x1 - Xi) 

i=l ' {i,j}eE(G) i=l 



We do so by the algorithm of Figure [5j In Figure [5] and in what follows, AB and C refer to parts of the 
equation as marked above. 

It is clear that the algorithm terminates for any finite graph, since the number of iterations of each for 
loop is finite. To prove correctness of the algorithm, we demonstrate that the following statement is true for 
each iteration of the for loop beginning in line 3: 

Prior to the m-th iteration, AB + C — 1 only contains terms of degree m + 1 or greater. 

(*) 

Furthermore, all terms of degree m + 1 are square-free. 

When the for loop is initialized, A is set to the polynomial — X)i=cP ^? ^(i,G) , and Qi,Qij are set to 
zero. Therefore, prior to the 0-th iteration, C = and AB + C — 1 = AB — 1 . Since the constant term in 
AB is equal to 

-CCnO,G)(-(a(G) + r))=~ a{G 1 ) + r (-(a(G)+r^j =1 , 

this implies that AB — 1 only contains terms of degree 1 or higher. Furthermore, linear terms are by definition 
square-free. Thus, (*) is true at initialization. 

Now we must show that if (*) is true prior to the m-th iteration of the for loop, then it will be true prior 
to the (m + l)-th iteration. Assume then, that prior to the m-th iteration, AB + C — 1 only contains terms 
of degree m + 1 or greater, and that all terms of degree m + 1 are square-free. We must show that prior to 
the (to + l)-th iteration, AB + C — 1 only contains terms of degree to + 2 or greater, and furthermore, that 
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all terms of degree m + 2 arc square-free. Prior to the m-th iteration, there are only two kinds of terms of 
degree to + 1 in AB + C — 1 : (1) terms corresponding to square-free, stable sets, and (2) terms corresponding 
to square-free, non-stable sets. We will show that both kinds of terms cancel during the m-th iteration. 

• Let x^Xfa ■ ■ ■x dm+1 be any square-free, stable set monomial in AB + C — 1 of degree (to + 1) . Since 
{di 7 d 2 , . . . ,d m+ i} is a stable set, all subsets of size to are likewise stable sets and appear as summands 
in ^(to, G) . Consider the coefficient of x dl x d2 ■ ■ ■ x dm+1 in C . During the m-th iteration, we apply rule 
1 and create this monomial ( m J ^ 1 ) times in C . Since this monomial is created by the multiplication of 
Xd 1 Xd 2 ---Xd^--- x dm+l with (x 2 dk -x dk ) , the coefficient for x dl x d2 ■ ■ ■ x dm+1 in C is 

1 Vm + i = -( TO + 1 ) C m + i 



Now we will calculate the coefficient of this same monomial in AB . This monomial is created in two 
ways, (1) multiplying -C^x dl x d2 ■ ■ ■ Xd j _ 1 £d j Xd j+1 ■ ■ ■ x dm+1 by x d] (repeated (m + 1) times, once for each 
x dj ), or (2) multiplying — C^ +1 x dl x d2 ■ ■ ■ x dm+1 by — (a(G) + r) (occurring exactly once). Therefore, the 
coefficient of x dl x d2 ■ ■ ■ x dm+1 in AB is 

-C% +1 {-(a(G)+r))-(m + l)Cg 

= a(G)C^ +1 + rC° +1 _^J^±^_^( a (G) + r - (m + 1)) 

= a(G)C° +1 + rCg +1 Cg +1 (a(G) + r (to + 1)) 

= a(G)C° +1 + rC° +1 - a(G)C% +1 rC G m+1 + (to + 1)C£ +1 

= (m+l)C£ +1 . 

Therefore, the coefficient for any square-free stable-set monomial in AB + C — 1 is 

(m+l)Cg +1 -(m + l)Cg +1 =0 . 

v " S v " 

from AB from C 

• Now, consider any square-free non-stable-set monomial x dl x d2 ■ ■ ■ x dm+1 in AB + C — 1 of degree to + 1 . 
Consider all ( m T ^ 1 ) subsets of {d\, di, ■ ■ ■ , d m+ i} , and let M be the number of stable sets among those 
("m X ) subsets. Each of those M subsets appears as a summand in Z?*{rn, G) . Therefore, the monomial 
x dl x d . 2 ■ ■ ■ x dm+1 is created M times in AB , and M times in C , by M applications of rule 2 . Therefore, 
the coefficient for any square- free non-stable set monomial x dl x d2 ■ ■ ■ x dm+1 in AB + C — 1 is 

from AB from C 

Finally, consider any non-square- free monomial x dl x d2 ■ ■ ■ x d ■ ■ ■ x dm+1 in AB + C — 1 of degree to + 2 . 
We note that {d\, ■ ■ . , d m+ i} is an stable set. To see this, note that every non-squarc-frcc monomial 
in AB is created by the product of an stable set with a linear term, and every non-square-free monomial 
in C is created by the product of an stable set with {x d — x dl ) for some di . During the m-th iteration, 
during applications of rule 1, x dl x d2 ■ ■ ■ x~d l ■ ■ ■ x dm+1 is added to Q dl . Therefore, when Q dl is subsequently 
multiplied by (x% — x dl ) , the monomial x dl x d ,, ■ ■ • x\ ■ ■ ■ %d m+1 is created. To summarize, this monomial is 
created in only one way in AB and only one way in C . Therefore, the coefficient for x^x^ • • • x^ • • • 
in AB + C - 1 is 

~Cm+l + pm+l = • 
from AB from C 

Therefore, we have proven that (*) is valid prior to every iteration of the for loop. Finally, upon termina- 
tion, i = a(G) + 1 . Therefore, by (*) , we know that AB + C — 1 only contains monomials of degree a(G) + 2 
or greater. However, there are no terms of degree a(G) + 2 in AB since deg( AB) = a(G) + 1 . Additionally, 
during the a(G)-th iteration, £?(a(G),G)xk is never an stable set. Therefore, only applications of rule 2 
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occur during the last iteration, and the final degrees of Qi, Qij are less than or equal to a(G) — 1 . Therefore, 
the monomial in C of greatest degree is of degree a(G) + 1 . Thus, there are no monomials in AB + C — 1 
of degree a(G) + 2 or greater, and upon termination, AB + C — 1 = . 
We have then shown that we can construct A, Qi, Qij such that 

-Y J C{i,G)^{i,G)\{~{a{G) + r)+Y. x A + Y Qa x i x i + Y & (*? " x <) • 

i=0 ' ^ i=l ' {i,j}e£(G) i=l 

* v ' V V ' v , 

A B c 

Since deg(Qj), deg(Q^) < (a(G) — 1) , and deg(A) = a(G) , this concludes our proof. □ 

Example 3.6. We display a certificate found by our algorithm. In Figure [6] is the Turan graph T(5,3) . 
It is clear that a(T(5, 3)) = 2 . Therefore, we "test" for a stable set of size 3. 

3 4 




Figure 6. Turan graph T(5, 3) 
The certificate constructed by our algorithm is 

1 = f - ~ (x 1 x 2 + x 3 x 4 ) - ~ (xi + x 2 + x 3 + X 4 + X 

1 1 

3 X2+ 3. 

n / j \ /I n /' i \ / 1 1 1 2 



1 


1\ 






3 , 2 H 


- g Jxix 3 + 




-i) 



j (xi + x 2 - 




CC4 - 


h x 5 - 3)+ 










- J 2:1X5 + 


Q x4 







jX 2 X4, + \^-jx 2 X 5 + l^-Xi + -JX 3 X 5 + \^-jX4X 5 + I -X 2 + -J (X 1 ~ Xi) + 

^1 + (-2 - x 2 ) + Qx 4 + (x§ - x 3 ) + Qx 3 + (2:4 - Xi) + f (a;§ - x 5 ) . 

Note that the coefficient for the stable set polynomial contains every stable set, and further note that 
every monomial in every coefficient is also indeed a stable set. 

We will now prove that the stability number a(G) is the minimum-degree for any Nullstcllcnsatz certificate 
for the non-existence of a stable set of size greater than a(G) . To prove this, we rely on two propositions. 

Proposition 3.7. Given a graph G , let M = {d\, d 2 , . . . , d|M|} & e an U maximal stable set in G . Let 



1 = A( - (a(G) + r) + ^2xi) + Qij^iXj + ^Q l {x 2 i - Xl ) 

^ «'=1 ' {i,j}£E(G) i=l 

be a Nullstellensatz certificate for the non-existence of a stable set of size a(G) + r (with r > 1), and let 

/ n \ n 

1 = A'( -( a (G)+r) + 5>i)+ Y Q'ioXiXj + Y.Q'ii.^-Xi) (3.6) 

, »'=1 , {i,0}eE(G) i=l 



c 



D 



be the reduced certificate via Lemma 3.4 Then, for i € {1, . . . , \M\\ , the linear term x^ appears in A' with 
a non-zero coefficient. 
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Proof. Our proof is by contradiction. Assume that x di does not appear in A' with a non-zero coefficient. 



By inspection of Eq. 3.6 we see that A' must contain the constant term ao = — (a(G) + r) . Therefore, the 
term aoX di appears in A' B . However, aoX di does not cancel within A'B since x di does not appear in A' by 
assumption. Therefore, aoX di must cancel with a term elsewhere in the certificate. Specifically, since C only 
contains terms multiplying edge monomials, a^x^ must cancel with a term in D . But the linear term aox^ 
is generated in D in only one way: ao must multiply {x 2 d . — Xd t ) Therefore, 

Q' d . = other terms + a . 

But when a multiplies (x d . — x di ) , this not only generates — agx^ , which neatly cancels its counterpart in 
A'B , but it also generates the cross-term aoX d . , which must cancel elsewhere in the certificate. However, 
x \, does not cancel with a term C , since C contains only terms multiplying edge monomials, and x 2 d . does 
not cancel with a term in A'B , since A' contains only terms corresponding to square-free stable sets and 
also because A' does not contain Xd i by assumption. Therefore, a,QX d . must cancel elsewhere in D . There is 
only one way to generate a second x d . term in D: agXdi must multiply (x d . — Xd t ) ■ Therefore, 

Q' d = other terms + OQXd\ + oq . 

Now, we assume 

Q'di = other terms + a,QX d . + aoa;^ -1 • ■ ■ + a x d . + agx^ + a . 

When aox d . multiplies (x d . — x^) , this generates a cross term of the form aox d +2 . This term must cancel 
elsewhere in the certificate. As before, aox^ +2 does not cancel with a term in C , since C contains only terms 
multiplying edge monomials, and aox^ +2 does not cancel with a term in A' B , since A' contains only terms 
corresponding to square-free stable sets. Therefore, aox d +2 must cancel elsewhere in D . But as before, there 
is only one way to generate a second aoX^ +2 term in D: ao£^ +1 must multiply {x d . — x di ) ■ Therefore, 

Q' d . = other terms + ao2^ +1 + aox d . + • • • + a^x 2 d . + a^Xd t + ao- 

To summarize, we have inductively shown that in order to cancel lower-order terms, we are forced to generate 
terms of higher and higher degree. In other words, Qd; contains an infinite chain of monomials increasing 
in degree. Since deg(Q^ ) is finite, this is clearly a contradiction. 

Therefore, Xdi must appear in A' with a non-zero coefficient. □ 

Proposition 3.8. Given a graph G , let M — {difd®, . . . ,d|M|} be any maximal stable set in G , and let 
{ci, C2, • • ■ , Cfc+i} be any (fc + 1) -subset of M with k < |Af| . Let 

(n \ n 

-(a(G) + r) + J2xA + E QijXiXj +^2Qi(xj - Xi) 
i=l ' {i,j}£E(G) t=l 

be a Nullstellensatz certificate for the non-existence of a stable set of size a(G) + r (with r > X), and let 

(n \ n 

-(a(G)+r) + J2xi]+ E Q-jXtXj + ^Qiix 2 - Xl ) 
i=l ' {i,j}eE(G) 1=1 
B c D 

be the reduced certificate via Lemma 3.4 If x Cl x C2 ■ ■ ■ x Ck appears in A' with a non-zero coefficient, then 



x Cl x C2 ' ' 'X Ck x Ck+1 also appears in A' with a non-zero coefficient. 

Proof. Our proof is by contradiction. Assume x Cl x C2 ■ ■ -x Ck with k < \M\ appears in A' with a non-zero 
coefficient, but x Cl x C2 ■ ■ -x Ck x Ck+1 does not. Since x Cl x C2 ■ ■ -x Ck appears in A' , x Cl x C2 ■ ■ -x Ck x Ck+1 clearly 
appears in A'B and must cancel elsewhere in the certificate. However, x Cl x C2 ■ ■ ■ x Ck x Ck+1 does not cancel 
with a term in C , since C contains only terms multiplying edge monomials and {ci, C2, ■ ■ ■ , c^+i} is a 
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stable set. Furthermore, x Cl x C2 ■ ■ -x Ck x Ck+1 does not cancel with a term in A' B , since A! does not contain 
x Cl x C2 ■ ■ ■ x Ck x Ck+1 by assumption. Therefore, x Cl x C2 • • • x Ck x Ck+1 must cancel with a term in D , and for at 
least one i , x Cl x C2 ■ ■ ■ x^ i ■ ■ ■ x Ck x Ck+1 appears in Q c . with a non-zero coefficient. 

When x Cl ■ ■ ■ x^ ■ ■ -x Ck x Ck+1 multiplies {x 2 . — x Ci ) , this generates a cross term of the form x Cl ■ ■ ■ x 2 . ■ ■ -x Ck x Ck 
which must cancel elsewhere in the certificate. Let mi = x Cl x C2 ■ ■ ■ x Ck x Ck+1 and mi — x Cl x C2 ■ ■ ■ x c . ■ ■ ■ x Ck x Ck+1 
Note that deg(m2) = deg(mi) + 1 . As before, mi does not cancel with a term in C , since C contains only 
terms multiplying edge monomials, and mi does not cancel with a term in A' B , since A' contains only 
terms corresponding to square-free stable sets. Therefore, mi must cancel elsewhere in D . 

In order to cancel mi in D , for some Cj , we must subtract one from the Cj-th exponent, and then multiply 
this monomial by [x 2 . — x Cj ) . However, this generates a cross-term UI3 = x Cl x C2 ■ ■ ■ x 2 . ■ ■ ■ x c ■ ■ ■ x Ck x Ck+1 
where deg(m3) = deg(m2) + 1 . Note that in the case when j = i , 7713 = x Cl x C2 ■ ■ ■ x^.. ■ ■ ■x Ck x Ck+1 , but 
deg(m3) still is equal to deg(m 2 ) + 1 . 

Inductively, consider the n-th clement in this chain, and assume it appears with a non-zero coefficient in 
some Q Cj . Let 

m — r ai T a2 ■ ■ ■ -r a k+l 
llt-n — X Cl X C2 X Ck + 1 , 

where a,- L > 1 for i G {1, . . . ,k+ 1} . When m„ multiplies (x\. — x Cj ) , this generates the cross-term m n x 2 . . 
This term must cancel elsewhere in the certificate. As before, m n x 2 . does not cancel with a term C , since 
C contains only terms multiplying edge monomials, and m n x\_ does not cancel with a term in A'B , since 
A' contains only terms corresponding to square-free stable sets and x Cl x C2 ■ ■ -x Ck x Ck+1 does not appear in 
A' by assumption. Therefore, m n x 2 . must cancel elsewhere in D . 
In order to cancel m n x 2 . in D , note that 

m r - 2 — T ai t Q2 • • • T a i+ 2 . . . rfO-k + l 

and for some I , let 

m , , — T "i r Q 2 . . . T aj+2 _ ai-l , , , a k + 1 
'«n+l — x Cl x C2 x c . x Cl ^ch+i ■ 

Note that deg(r7i Jl+ i) = dcg(rn„) + 1 . Therefore, in order to cancel m n x 2 . , m n +i multiplies (x 2 t — x Cl ) , 
which generates a new term of higher degree: m n +\x\ . 

To summarize, we have inductively shown that in order to cancel lower-order terms, we are forced to 
generate terms of higher and higher degree. In other words, mi, m%, . . . form an infinite chain of monomials 
increasing in degree. Since deg(Q9 is finite, this is clearly a contradiction. 

Therefore, x Cl x C2 ■ ■ ■ x Ck x Ck+1 must appear in A' with a non-zero coefficient. D 

Using Propositions |3.7| and [378] we can now prove the main theorem of this section. 



Theorem 3.9. Given a graph G , any Nullstellensatz certificate for the non-existence of a stable set of size 
greater than a(G) has degree at least a{G) . 

Proof. Our proof is by contradiction. Let 

(n \ n 

-(a(G) + r) + J2 x i) + E QijXiXj + J2Q i (x 2 i - Xi) 
i=l ' {i,j}€E(G) i=l 

be any Nullstellensatz certificate for the non-existence of a stable set of size a(G) + r, with r > 1 , such that 
deg(A),deg(Q i ),deg(Q ij ) < a(G) , and let 

(n \ n 
-(a(G)+r) + ^2x i )+ Yl Q',, ■Y.Q'^: •<•-) (3-7) 
1=1 ' U,j}eE(G) i=l 
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be the reduced certificate via Lemma 3.4 The proof of Lemma 3.4 implies deg(A') < deg(A) < a(G) . Let 
M = {di, d 2 , • ■ ■ 7 ^q(g)} be any maximum stable set in G . Via Proposition 3.7 we know that x^ appears in 



A' with a non-zero coefficient, which implies (via Proposition 3.8 1 that x^x^ appears in A' with a non-zero 
coefficient, which implies that x^x^ Xd 3 appears in A' and so on. In particular, x^x^ ' ' ' x d a i G ) appears 
in A' . This contradicts our assumption that deg(A') < a(G) . Therefore, there can be no Nullstcllensatz 
certificate with deg(A) < a(G) , and the degree of any Nullstellensatz certificate is at least a(G) . □ 

Propositions |3.7| and |3.8| also give rise to the following corollary. 



Corollary 3.10. Given a graph G , any Nullstellensatz certificate for the non-existence of a stable set of 
size greater than a(G) contains at least one monomial for every stable set in G . 



Proof. Given any Nullstellensatz certificate, we can create the reduced certificate via Lemma 3.4 The 



proof of the Lemma 3.4 implies that the number of terms in A is equal to the number of terms in A 1 . Via 



Propositions 3.7 and 3.8 A 1 contains one monomial for every stable set in G . Therefore, A also contains 



one monomial for every stable set in G . D 

This brings us to the last theorem of this section. 

Theorem 3.11. Given a graph G , a minimum- degree Nullstellensatz certificate for the non-existence of a 
stable set of size greater than a(G) has degree equal to a(G) and contains at least one term for every stable 
set in G . 



Proof. This theorem follows directly from Theorems 3.5 3.9 and Corollary 3.10 □ 



Finally, our results establish new lower bounds for the degree and number of terms of Nullstcllensatz 
certificates. In earlier work, researchers in logic and complexity showed both logarithmic and linear growth 
of degree of the Nullstcllensatz over finite fields or for special instances, e.g. Nullstellensatz related to the 
pigeonhole principle (see [5] , and references therein) . Our main complexity result below settles a question 
of Lovasz 24J : 

Corollary 3.12. Given any infinite family of graphs G n , on n vertices, the degree of a minimum- degree 
Nullstellensatz certificate for the non-existence of a stable set of size greater than a(G) grows as f2(n) . 
Moreover, there are graphs for which the degree of the Nullstellensatz certificate grows linearly in n and, 
at the same time, the number of terms in the coefficient polynomials of the Nullstellensatz certificate is 
exponential in n . 

Proof. The stability number of a graph G with n nodes and m vertices grows linearly ([H]) since 

a(G) > ^ ^(2m + n+ 1) - ^(2m + n+l) 2 - 4n 2 ^j . 

Finally, it is enough to remark that there exist families of graphs with linear growth in the minimum degree 
of their Nullstellensatz certificates, but exponential growth in their numbers of terms. The disjoint union 
of n/3 triangles has exactly 3™/ 3 maximal stable sets. Therefore, its Nullstellensatz certificate's minimum 
degree grows as 0(n/3) , but its number of terms grows as 2> n / 3 (see [T3] and references therein). □ 
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3.2. The Nullstellensatz and 3-colorability 

In this section, we investigate the degree growth of Nullstellensatz certificates for the non-3-colorability for 
certain graphs. Curiously, every non-3-colorable graph that we have investigated so far has a minimum-degree 
Nullstellensatz certificate of degree four. Next, we prove that four is indeed a lower bound on the degree of 
such certificates. 



Theorem 3.13. Every Nullstellensatz certificate of a non- 3- colorable graph has degree at least four. 



Proof. Our proof is by contradiction. Suppose there exists a Nullstellensatz certificate of degree three or 
less. Such a certificate has the following form 



1 = E P W ( x * ~ *) + E P {y > + x * x i + x ^ ' 



(3.8) 



1=1 



where P{j} and Puj} represent general polynomials of degree less than or equal to three. To be precise, 

n n n 

P M = E Hi}^ + E E b {i}*t X l X t 



S = l t=l 



X! X! X! C{i} stu X s X t X u + ^ E d {i}stX s X t + ^ e {i}s x s + f{i} 



s=l t=s+l u=t+l 



8 = 1 t=l 



and 



p m=J2 a 



wis 



EE 6 

s = l t=l 



{ij}stX s Xf 



+ C{.y} stu a; s a;ta; n + d^j^ st x s Xt + e{jj-} s a; s + /{y} ■ 

s=l t=s+l «=t+l s=l t=l s=l 

Since we work with undirected graphs, note that a{ij} s = o-{ji} s > an d this fact applies to all coefficients a 
through / . Note also that when is not an edge of the graph, P^ — and thus a,u,j\ s = . Again, this 

fact holds for all coefficients a through / . 

When P{jj multiplies (xf — 1) , this generates cross-terms of the form and —P{i} ■ In particular, this 

generates monomials of degree six or less. Notice that P^j^(x 2 + 2 ) does not generate monomials of 

degree six, only monomials of degree five or less. We begin the process of deriving a contradiction from Eq. 
by considering all monomials of the form x^xf that appear in the expanded Nullstellensatz certificate. 



3.8 



These monomials are formed in only two ways: Either (1) x^(xf — 1) , or (2) xf(x^ — 1) . Therefore, the 
equations for x^xf are as follows: 



= 
= 



a {l}l ! 

0{1}2 + a {2}l , 



(coefficient for x\x\ = xf) 
(coefficient for x\x\) 



(7.1) 
(1.2) 



= 
= 



fl{n-i}n + a{„}(„_i) , (coefficient for xf n _^x^) 
fl{n}n ■ (coefficient for x n x n = x n ) 



Summing equations 1.1 through I.n 2 , we get 



n n 

= EE a ws 

1 = 1 s=l 



(I.n 2 - 1) 
(I.n 2 ) 

(3.9) 
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Let us now consider monomials of the form x 2 s x t x\ (with s / i). These monomials are formed in only one 
way: by multiplying b^ st x 2 x t by x\ . Therefore, since the coefficient for x 2 s x t x\ must simplify to zero in 
the expanded Nullstellensatz certificate, bu\ st — for all b^ . When we consider monomials of the form 
x s x t x u xf (with s < t < it) , we see that c^ stu — for all , for the same reasons as above. 

As we continue toward our contradiction, we now consider monomials of degree three in the expanded 
Nullstellensatz certificate. In particular, we consider the coefficient for x^ . The monomial x^. is generated in 
three ways: (1) f{ s }(x^. — 1) , (2) au\ s Xg (x? — 1) (from the vertex polynomials), and (3) e^ st y s x s (x 2 +x s Xt+x 2 ) 



(from the edge polynomials). The equations for a^, 



c~ are as follows: 



= 



= 



f{i} -^2 a { l }i+ E e {i*}i ' (coefficient for x J) (II. 1) 

i=l teAdj(l) 



/{2} - E a W 2 + E e { 2t > 2 ' (coefficient for x%) 

i=l teAdj(2) 



(17.2) 



n 



«= /w-E 

Summing equations 77.1 through 77. n , we get 



a {i}n + E e {™*}^ 

= 1 iGAdj(n) 



(coefficient for a^) 



(77.n) 



= E /{<> - ( E E a w s ) + E E e ( s *} S 

i=l ^ 1=1 s = l ' s = l teAdj(s) 



n n 



(3.10) 



Since the degree three or less Nullstellensatz certificate (Eq. 3.8) is identically one, the constant terms must 
sum to one. Therefore, we know X)"=i /{«} = — 1 ■ Furthermore, recall that e^ s t\ s — if the undirected edge 
{s,t} does not exist in the graph. Therefore, applying Eq. 3.9 to Eq. 3.10 we have the following equation 



1 - E e l st l s 

s,t=t, 



(3.11) 



To give a preview of our overall proof strategy, the equations to come will ultimately show that the right-hand 



side of Eq. 3.11 also equals zero, which is a contradiction. 

Now we will consider the monomial x 2 Xt (with s ^ t). We recall that b^ st = for all b^y (where b^y st 
is the coefficient for x 2 Xt in the i-th vertex polynomial). Therefore, we do not need to consider b^ st in the 
equation for the coefficient of monomial x 2 x t . In other words, we only need to consider the edge polynomials, 
which can generate this monomial in two ways: (1) e^ st } s x s ■ x s Xt , and (2) e^ s ^ t Xt ■ x 2 . The N = 2(") 
equations for these coefficients are: 



(coefficient for x\x<i) 



e{i2}i + E e {i«}2 ■ 
ieAdj(i) 

e {i3}i + E! e {i«}3 ' (coefficient for x\x%) 
ieAdj(i) 



(777.1) 
(777.2) 



e{n(n-l)}n + 



E 

iSAdj(r, 



'{ni} [n— 1) 



(coefficient for x 2 n x n _\) 



(777.7V) 
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When we sum equations 1 1 1. 1 to 1 1 I.N , we obtain 



EE e {^ + (E E e ^}*) + (E E E e ^}"H°- ( 3 - 12 ) 



=1 t=l, x s=l t^Adj(s) 7 v s=l teAdj(s) u=l, 



partial sum A y 

partial sum B 

However, recall that e^ st y u — when {s,t} does not exist in the graph, and also that e^ st y t — es ts x t . Thus, 
we can rewrite partial sum A from Eq. |3.12| as 

n n n n n 

E E e {st}* = EE e { s *}« = EE e { ts }* ■ 

s=l teAdj(s) s=l t=i, s=l t=i. 

Substituting the above into Eq. 3. 12| yields 



2 E e {^ + (E E E e ^}u)=o- (3-i3) 

.,*=i, v s=iteAdj(s) "=i. y 

v ' 

partial sum B 

Finally, we consider the monomial x s xtx u (with s < t < u). We have already argued that c^ stu — for all 
C{ij (where cuy s t u is the coefficient for x s xtx u in the z-th vertex polynomial). Therefore, as before, we need 
only consider the edge polynomials, which can generate this monomial in three ways: (1) es st \ u x u • x s Xt , (2) 
e{su}tXt ■ x s x u , and (3) est u \ s x s ■ Xtx u . As before, these coefficients must cancel in the expanded certificate, 
which yields the following (™) equations: 

0= e {12}3 + e {13}2 + e {23}1 , (coefficient for XiX 2 x 3 ) (IV.l) 

= e{i2}4 + 6{i4}2 + e{24}i , (coefficient for X1X2X4) (-^-2) 



e{(«-2)(n-i)}n + e{( n _ 2 )n}(n-i) + e{(„-i)n}(«-2) • (coefficient for x n _ 2 x„_ix„) (IV.M) 



Summing equations IV. \ through IV.M , we obtain 

n— 2 n— 1 H / \ 

E E E ( e {st}u + e{su}t + e{tu}s j = ■ (3.14) 

s=l t=s+l «=t+l ^ ' 

Now we come to the critical argument of the proof. We claim that the following equation holds: 

1 n \ / n— 2 Ti—l n / \\ 

E E E e M» = 2 E E E l e {st}u + e {su}t + e {tu}s ) ) . (3.15) 



s=l te_Ad/(s) 11=1, 



s=l t=s + l u=t+l 



Notice that the left-hand and right-hand sides of this equation consist only of coefficients e^ st j u with s,t,u 
distinct. Consider any such coefficient &s s t\u . Notice that e/ st T. u appears exactly once on the right side of the 
equation. Furthermore, either e^ st j u appears exactly twice on the left side of this equation (since s € Adj(i) 
implies t € Adj(s)), or e{ st } u = (since the edge {s,t} does not exist in the graph). Therefore, Eq. 3.15 is 
proven. Applying this result (and Eq. 3.14) to Eq. 3.13|gives us the following: 



E H 



st}s 



l<»,t<r. 



But Eq. 3.16 contradicts Eq. 3.11 (1 = 0), thus there can be no certificate of degree less than four. 



(3.16) 



□ 
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It is important to note that when we try to construct certificates of degree four or greater, the equations 
for the degree-6 monomials become considerably more complicated. In this case, the edge polynomials do 
contribute monomials of degree six, which causes the above argument to break. 

We conclude this subsection with a result that allows us to bound the degree of a minimum-degree 
Nullstellensatz certificate of a particular graph, if that graph can be "reduced" to another graph whose 
minimum-degree Nullstellensatz certificate is known. 

Lemma 3.14. 

1. If H is a subgraph of G , and H has a minimum- degree non- 3 -color ability Nullstellensatz certificate of 
degree k , then G also has a minimum-degree non-3-colorability Nullstellensatz certificate of degree k . 

2. Suppose that a non- 3 -colorable graph G can be transformed to a non- 3 -colorable graph H via a sequence 
of identifications of non-adjacent nodes of G . If a minimum- degree non- 3 -color ability Nullstellensatz 
certificate for H has degree k , then a minimum- degree non- 3 -colorable Nullstellensatz certificate for G 
has degree at least k . 

Proof. 

Proof of 1: Since H is a subgraph of G , then any Nullstellensatz certificate for non-3-colorability of H is 
also a Nullstellensatz certificate for non-3-colorability of G . 

Proof of 2: Suppose that G has a Nullstellensatz certificate for non-3-colorability of degree less than k . 
The certificate has the form 1 = ^a^u,; + X)^{u'} e {ij} w here = xf — 1 , — xf + XiXj + Xi , and 

both ai and b^jj denote polynomials of degree less than k . Since the certificate is an identity, the identity 
must hold for all values of the variables. In particular, it must hold for every variable substitution Xi = Xj 
when the nodes are non-adjacent. In this case, the variable reassignment (pictorially represented in Figure [7]) 
yields a Nullstellensatz certificate of degree less than k for the transformed graph H . Note that the parallel 
edges that may arise are irrelevant to our considerations (see such examples in Figure^). But this is in 




Figure 7. Converting G (the 5-odd-wheel) to H (the 3-odd-wheel) via node identifications. 

contradiction with the assumed degree of a minimum-degree certificate for H . Therefore, any certificate for 
G must have degree at least k . □ 

3.2.1. Cliques, Odd Wheels, and Nullstellensatz Certificates 

Theorem 3.15. For K n with n > 4 , a minimum- degree Nullstellensatz certificate for non-3-colorability 
has degree exactly four. 

Proof. It is easy to see that K 4 is a subgraph of K 5 , which is a subgraph of K e , and so on. The decision 
problem of whether K4 is 3-colorable can be encoded by the system of equations 

x\ — 1 = 0, X2 — 1=0, x\ + x\x 2 + x\ — , x\ + £1X3 + x\ = , x\ + X1X4 + x\ = , 
x\ — 1 = 0, x\ — 1=0, x% + X2X3 + x§ = , x\ + x 2 X4 + x\ = , x\ + x 3 Xi + x\ = , 
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with one equation per vertex and one equation per edge. Using the linear-algebra heuristic described at the 
beginning of Section [3] we find that K4 has a degree-4 Nullstellensatz certificate for non-3-colorablility: 



(— x\ — 1)(CE? - 1) + (^ x 4 ~ ^ X 4 X 2 - ^ X 4 X 3 - ^ X 4 X 1 + ^X%X2Xl + ^X^X^X^j 



X%X\ + -x 4 x 2 xi + -x^x^x-ij (x 4 + X2X4 + x 2 ) 

-x\ + \x\x2 - -x\x\ - \x\x2xx J (x\ + x 3 x 2 + x\) + \ \x\ + \x\x 2 + \x%x\ + \x\x%x\ J {x\ + x 3 x 4 + xf) 
9 9 9 9 / \9 9 9 9 / 

->- x\x\ - x 4 x\ + x\\ (x\ + x\X4 + x\) + \x\x 2 {x\ + x\X2 + x\ ) + ( - \x\ - \x\x2 \ ix\ + X\X 3 + x\ ). 



3 4 ■ , - " • 3 • - ' ' ' - ' V 3 - 3 



Because K4 has a degree-4 Nullstellensatz certificate as shown above, K n for n > 4 also has a degree-4 



Nullstellensatz certificate via Lemma 3.14 (1). □ 



The odd- wheels consist of an odd-cycle rim, with a 
rather easy to see that no odd- wheel is 3-colorable. It 
degree Nullstellensatz certificate for non-3-colorablility. 



center vertex connected to all other vertices. It is 
is natural to ask about the degree of a minimum- 



Theorem 3.16. Given any odd-wheel, the degree of a minimum- degree Nullstellensatz certificate for non- 
3-colorability is four. 



Proof. First, we will prove that there exists a certificate of degree four for the n-th odd- wheel. We begin 
by displaying a degree-4 certificate for the 3-odd-wheel: 



1 = 



~x\ - ^x\x 2 - ^x\x 3 - ^x\x + ^xlx 2 x + ^xlx 3 xo ) {x{ + X 2 Xi + x 2 )+ 



1 A , 2 J, 



1 



9 



9 



9 



1 



x 1 + -x x x 2 - -x\x - -x^Xq ) (x 2 + x 3 x 2 + x 3 ) + -x 1 x 2 (x 2 + x Q x 2 + x ) + 



9 



9 



x 1 + -x Y x 2 + -x^q + -x^xq ) (x l + x 3 xi + x 3 ) + -x 1 {x 1 + x xi + x )+ 



:xi 



1 



9 



\x 2 (x§ + x x 3 + x z Q ) + {-x\ - l)(xf - 1). 



(3.17) 



3 1 3 

For now, we denote the non-3-colorability certificate for the 3-odd-wheel as follows: 

1 = "{12} e {12} + a {23} e {23} + "{20} e {20} + «{13} e {13} + "{10} e {10} + «{30} e {30} + "l^l, 

where V\ = x\ — 1 , and e^jj = x\ + XiXj + x 2 and a\ and ot^j} denote polynomials of degree four in 
M.[xq, Xi, x 2l x 3 ] . In Figure [8j we can see that the topological difference between the 3-odd-wheel and the 




Figure 8. Here we show the evolution of the 3-odd-wheel to the 5-odd-wheel. This can be extended from the 

n-th odd wheel to the (n + 2)-th odd wheel. 

5-odd-wheel is that the edge (1, 3) is lost, and the vertices 4, 5 and associated edges (3, 4), (4, 5), (5, 1), (0, 4) 
and (0, 5) are gained. We can exhibit an algebraic relation (or syzygy) as follows: 

a {13} e {13} = «{13} e {15} + /3{34>e{34} + /?{45} e {45} + /3{01} e {01} + P{03} e {03} + /3{04>e{04} 
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/3{05} e {05} j 



(3.18) 



where G ^\xq, Xi, X2, £3, X5] and deg(/3/jj}) < 4 . Note that the coefficients for e{i3j and e{i5j are 

the same. Recall that 



«{13} 



1 



^x\x + ^xfx 2 x 



Since a{ 13 } does not contain the variable X3 , «{i3} € M[xq, X\, x%\ and not in M[xo, 11, X2, £3] ■ Therefore, 
because there exists a degree-4 certificate for the 3-odd-wheel, we can simply use the above syzygy (Eq. 



3. 18 1 to substitute for the a{ 13 }e{ 13 } term in the degree-4 3-odd-whccl certificate (Eq. 3.171. The resulting 



polynomial is a degree-4 certificate for the 5-odd- wheel where the coefficient for e{ 15 j in the 5-odd- wheel 
certificate is exactly the same as the coefficient for e{ 13 } in the 3-odd-wheel certificate (both coefficients are 
a {i3})- Thus, we can again use the syzygy of Eq. 



3.18 



(with the variable substitutions of £4 



Xfy ; X§ 



x 7 



and X3 — ► X5), and substitute for the a{ 13 }e| 15 } term in the 5-odd-wheel certificate to obtain a degree-4 
7-odd-whcel certificate. Thus, by induction, we obtain degree-4 certificates for all odd-wheels. It remains for 
us to show that such a syzygy exists. 

The special syzygy was found via the linear algebra heuristic described at the beginning of Section |3]and is 
listed below. Note that the coefficients for e{ 13 j. and are indeed identical and are equal to the coefficient 
for e{ 13 } in the 3-odd-wheel certificate we presented earlier (Eq. 3.171. 



= - 



~x\ + \x\x2 + \x\x + ^xfx 2 x ) (x\ + X3X1 + x\) + \ \x\ + \x\x2 + \x\x + \x\x2X$ ) (x\ + X 5 Xl 
9 9 9 / v ^ ' \ 9 9 9 9 / •> ^ 

7 V ' e t 13 > " T" ' e f 15 > 



-4) 



"{13} 



^x?x + ixix 2 x X5 - \x\X2X A xc, - ^-x 1 x 3 Xq - ^-x 1 x 3 x X4 - -x^o - \x2x\x4 + \x\ ) {x\ + X3X4 + x|) 
v 9 9 9 9 9 99 9/ s ^ / 



"{13} 
1 



e {34} 



2 4 2 o 



L 9 1 9 1 

-XJX2X4 H XJX0X4 X1X2X3XQ ■ 

) 9 9 



1 1 2 , 1 2 2 4 

-XlX 2 X 3 X 4 - -x 1 x 2 x + -xix 2 x 4 - -x 



1 



1 



g X X4 — -x 4 



\x4X5 - \xaxI ) (x\ + X4X5 + xl) 
1 9 / > v s 



= {45} 



3 2 "^2 2 2 

"'o " . 1 — X\X^XQ — — X-j^X3X0 ~\~ — X-^X4Xfj 

9 9 3 9 



2 2 1 2 

— XQX4X5 — -X1X4X5 
2 



1 2 

-X3XOX4 



X2X 3 x\ 



^X\X 3 X% - ^X^XOxl 

- ^xfx 2 x 3 - ix?x 2 x 5 + ^xfx 3 - ^xfx 5 + ixfx x 5 - ?i 

5 O 4 rjl r) X n 2 Q 2 n l_o 1 Q 2 o 1 

— x\Xqx\ xxXqx^ x\x\x^ 3:i£q H — X2X3XQ H — X2X3X4 2:23:33:5 H — x^xqx\ H — a:2a:3a:o a:; 4 

9999999993 

1 



rrjrrp + -^x\x\ — -X\X^X4 — — X1X3X0X4 — — X\XqX4X§ 



1 1 o 

-0:2X3^03:5 + -X2X4 - 



x\xq 



1 4 1 3 

-zj 1; 



/2 4 1 3 4 3 4 3 1 
" I g^l + g^l^ + -x^o + -x ia ;4 - - 



9 



X4X2X4 



X 3 X4 



2 o 2 



T: x o x i ~ \ x o x a ) ( x o + x o x i + x\) 
9 9 / x . ' 



e {01} 



1 



2 2 



2 2 1 1 1 

-X1X2X Q — -X1X2XQX4 — -X1X2XQX5 + - 



2 
9 

1 3 1 2 Z 

-x 3 x - -x 3 x x 4 - - 



3 

X1X2X4X5 



L 2 

-X^XQ 



-X1X 3 X - 



9 



Xi X 3 X4 



9 



2 2 
■ x i x o " 



9 



x\xqX4 



2 2 2 



2 2 1 

-xia;3a; + -a;ia;3a;oa:4 ■ 



22 
Z X 3 X 



*-x 3 x x 2 4 - -x^ - -x 3 X4 ) (xl + x x 3 + xl) 
9 9 9 / v ^, ✓ 



£ {03} 



9 X 1 X 5 ~ ^ x l x 2 x 3 



1 



9' 



x i x 2 x 5 - - x i x l - 

^X 1 Xq + ^X 1 XqX4 + ^X 1 X( I X 5 + ^x 2 x 3 x x 5 



1 1 2 1 2 1 2 2 

-XlX 2 X 3 X4 + -XlX 2 X - -X1X2X4 + -XlX 3 X Q + -XlX 3 X X4 



1 



L 2 

-x 2 x 3 x 5 - 



2 3 
■ y x 3 x 



-xlx X4 



L 2 2 

-2:3X4 



1 3 

3Z3Z0 



-x 3 xl 



- -x 3 x x 4 

1 3 1 3 1 

-x{x 2 + -x^x 4 + - 



-Xq ) (Xq + X0X4 + xl) 
i J ' v ' 



e {04} 



x\x2X 3 



- X -x\x 2 X4 - ixfxg - 



1 1 

XlX 2 X 3 Xo XlX 2 X 3 X 4 + - 

1 9 9 



g xix 2 x - -xix 2 x 4 
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- -SiZq + -x-iXqX4 - -x 2 x 3 x X4 - -X2X3X4 - -x x\x 5 - -x x 4 xl + -x\x\ + -x±x\ ) (Xg + x x 5 + x\) . 
99 9 9 9 9 99/ \ v ✓ 

□ 

Finally, the reader may easily observe that Theorem |1.2[ Part 2 follows directly from Theorem |3.15| 
Theorem |3.16[ and Lemma [3.14| 

3.3. Nullstellensatz Certificates for Other Non-3-Colorable Graphs 

With the aid of a computer, we searched many families of non-3-colorable graphs, hoping to find explicit 
examples with growth in the certificate degree. Every graph we have investigated so far has a Nullstellensatz 
certificate of degree four. This suggests that examples with degree growth are rare for graphs with few 
vertices, and that many graphs have short proofs of non-3-colorability. In Figure [9] we describe the Jin and 
Grotzch graphs, and in Figure [lOj we describe the "Flower" family. Kneser graphs are described in most 
graph theory books. In Table [l] we present a sampling of the many graphs we tried during our computational 
experiments. Note that we often used our probabilistic linear algebra algorithm, selecting p = .4 as a likely 
threshold for feasibility 



Graph 


vertices 


edges 


row 


col 


V 


dog 


flower 8 


16 


32 


51819 


49516 


. 1 


4 


flower 10 


20 


40 


178571 


362705 


1 


4 


flower 11 


22 


44 


278737 


278844 


.5 


4 


flower 13 


26 


52 


629666 


495051 


.4 


4 


flower 14 


28 


56 


923580 


705536 


.4 


4 


flower 16 


32 


64 


1979584 


1674379 


.4 


4 


flower 17 


34 


68 


2719979 


2246535 


.4 


4 


flower 19 


38 


76 


4862753 


3850300 


.5 


4 


kneser-(6,2) 


15 


45 


39059 


68811 


.5 


4 


kneser- (7, 2) 


21 


105 


230861 


558484 


.5 


4 


kneser-(8,2) 


28 


210 


1107881 


3307971 


.5 


4 


kncser-(9,2) 


36 


378 


1107955 


3304966 


.5 


4 


kneser-(10,2) 


45 


630 


15,567,791 


36,785,283 


.5 


4 


jin graph 


12 


24 


12168 


13150 


.4 


4 


Grotzsch 


11 


20 


7903 


8109 


.4 


4 


G + {(3,4)} 


12 


24 


12,257 


13,091 


.4 


4 


G + {(7,12)} 


12 


24 


12,201 


13,085 


.4 


4 


G + {(1,8)} 


12 


24 


12,180 


13,124 


.4 


4 


G + {(3,4),(12,7)} 


12 


25 


12,286 


13,804 


.4 


4 



Table 1 .Experimental investigations for Flowers, Kneser, the Jin graph and the Grotzch graph. Here G 
denotes the uniquely-colorable graph displayed in Figure [9] 

A uniquely 3-colorable graph is a graph that can be colored with three colors in only one way, up to 
permutation of the color labels. Figure [9] displays a uniquely 3-colorable triangle-free graph [BJ. Since the 
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Figure 9. These graphs (from left to right) are (1) a uniquely 3-colorable graph, labeled with its unique 
3-coloring [6], (2) the Grotzsch graph, and (3) the Jin graph. 

graph is uniquely 3-colorable, the addition of a single edge between two similarly-colored vertices will result 
in a new non-3-colorable graph. Table [T] also details these experiments. Finally, we investigated all non-3- 
colorable graphs on six vertices or less: every one has a Nullstellensatz certificate of degree four. 




Figure 10.3,4 and 5 flowers (left to right). Note that the 3-flower is 3-colorable, whereas the 4 and 5 flowers 
are non-3-colorable. It is easy to see that only flowers that are multiples of 3 are 3-colorable. 
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